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The phase diagram of surface structures for the model catalyst RuO2(110) in contact with a gas 
environment of O2 and CO is calculated by density-functional theory and atomistic thermodynamics. 
Adsorption of the reactants is found to depend crucially on temperature and partial pressures in 
the gas phase. Assuming that a catalyst surface under steady-state operation conditions is close 
to a constrained thermodynamic equilibrium, we are able to rationalize a number of experimental 
findings on the CO oxidation over RuO2(110). We also calculated reaction pathways and energy 
barriers. Based on the various results the importance of phase coexistence conditions is emphasized 
as these will lead to an enhanced dynamics at the catalyst surface. Such conditions may actuate an 
additional, kinetically controlled reaction mechanism on RuO2(110). 
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I. INTRODUCTION 

A prerequisite for analyzing and understanding the 
electronic properties and the function of surfaces, is the 
detailed knowledge of the atomic structure, i.e., the sur- 
face composition and geometry. In this respect the ex- 
perimental techniques of ultra-high vacuum (UHV) Sur- 
face Science have significantly helped to build our cur- 
rent understanding^ and sometimes these UHV results 
can be related to high-pressure applications, like catal- 
ysis or corrosion. Often, however, such extrapolation is 
not straightforward and may even be impossible. 2 This is 
due to the fact, that some surface structures that can be 
stabilized in UHV - sometimes by sophisticated prepa- 
ration, annealing, and post-dosing procedures - may not 
exist at high temperature, T, and high pressure, p. And 
structures that exist at high pressure may be hard or 
even impossible to identify or stabilize under UHV con- 
ditions. Obviously, knowing what structures dominate 
the high-pressure application one is interested in, is cru- 
cial for finding a way to prepare the same or a similar 
situation under UHV conditions, in order to then per- 
form controlled, atomistic experiments. Thus, the deter- 
mination of a (T, p) phase diagram, covering the surface 
phases from UHV to realistic conditions appears to be 
critical for a meaningful and safe bridging of the pres- 
sure gap, enabling a Surface Science study relevant to 
the high-pressure problem. 

The concept of first-principles atomistic 
thermodynamic&££& enables us to calculate such 
surface phase diagrams, and for metal oxides this 
approach has recently proven to be most valuable (see 
for example Refs. l7l8l l MHy> . Using density- functional 
theory, one calculates the free energies of all plausible 
surface compositions and geometries, in order to identify 
the lowest-energy structure for a given condition of the 
thermodynamic reservoirs for the atoms and electrons. 
In Ref. we described the approach in detail for a 



one-component gas phase, and we now extend it to 
a multi-component environment also employing the 
concept of a "constrained thermodynamic equilibrium" ; 
at the end we will also analyze (and emphasize) the role 
of kinetic effects. A short paper on this study already 
appeared^ 

We believe that the below discussion is relevant to oxi- 
dation catalysis at a ruthenium catalyst, but the applica- 
bility of the methodology is much wider. From the calcu- 
lated surface phase diagram, a number of conclusions on 
the experimentally reported high efficiency of this oxide 
surfaco 11 ! 12 ^ 3 ! 14 ^ 5 ! 16 ^ 7 ! 18 ^ 9 can already be drawn. Most 
notably, we emphasize that gas phase conditions, that 
correspond to a coexistence of different surface phases, 
may lead to an enhanced dynamics and be particularly 
important for catalytic applications. For the specific ex- 
ample of the RuO2(110) surface we find that such phase 
coexistence conditions at high partial pressures enable a 
reaction mechanism, that does not play a role otherwise. 
Interestingly, the energy barrier is found to be even lower 
than that of other mechanisms considered so far«i£*22i 



II. THEORY 

Our approach connects density-functional the- 
ory (DFT) total-energy calculations and atomistic 
thermodynamics ^ ! i 4 i 5 : 6 In a preceding publication the 
procedure was described for an oxide surface in equilib- 
rium with a one-component gas phase 9 . We will therefore 
recapitulate it here only briefly, concentrating now on 
the extension to a multi-component gas phase and a 
"constrained thermodynamic equilibrium" . For clarity 
we will stick to the specific example of a RuO2(110) 
surface in contact with a gas phase environment formed 
of O2 and CO. The generalization to other surfaces in 
contact with arbitrary multi-component gas or even 
liquid phases is straightforward. 
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A. Surface free energy 

For a surface in equilibrium with atomic reservoirs (de- 
fined by a gas or liquid phase environment, or a macro- 
scopic bulk phase) , the most relevant structures are char- 
acterized by a low surface free energy, which is defined 

as 



7(T,fe}) = ^ 



G - ^2NiiM(T,pi 



(1) 



Here, G is the Gibbs free energy of the solid with the 
surface we like to study. If a slab is used, there are two 
surfaces (the top and the bottom side), both of which 
are then, of course, to be considered in the (total) sur- 
face area, A. jLtj(T,pj) is the chemical potential of the 
species of «th type (here: i = Ru, O2, CO), and Ni is the 
number of atoms (or molecules) of the ith species in the 
considered solid. T and {pf\ are the temperature and 
the partial pressures of the various species. Consider- 
ing the presence of two independent reservoirs of O2 and 
CO implies that O2 and CO molecules, though their mix 
forms the environment, are not in equilibrium with each 
other: From energy considerations alone, CO2 would re- 
sult as the most stable gas phase molecule for almost all 
temperature and pressure conditions, if the environment 
were able to attain thermodynamic equilibrium with it- 
self. However, the large free energy barrier for the gas 
phase reaction CO + 1/2 O2 — > CO 2 prevents that this 
reaction plays any role on the time scales of interest. 
Thus, our treatment ignores (because of good reasons) 
CO2 formation in the gas phase, and only at the end of 
our study (sections III D-F) we will consider that such 
a reaction may take place between species that are ad- 
sorbed on the surface. As a consequence, the oxygen and 
CO chemical potentials in this "constrained thermody- 
namic equilibrium" situation are given by the expressions 
discussed in the next paragraph. 

With respect to the Ru chemical potential we note the 
other constraint of the present study, namely the pres- 
ence of macroscopic quantities of bulk Ru02 . If the tem- 
perature is not too low, this Ru02 material is in equilib- 
rium with the O2 environment which implies 

MRu + MO2 = feuOji ( 2 ) 

where g R ^ l k (T, pr u o 2 ) is the Gibbs free energy of the 
bulk oxide (per formula unit). Inserting this into eq. Q 
to eliminate /zr u and rewriting it for a slab calculation 
with two equivalent surfaces leads to a surface free en- 
ergy, which is now a mere function of the chemical poten- 
tials determined by the equilibrium with the gas phase 
reservoirs, namely fico and /io 2 > 



l(T,Po 2 ,Pco) = -j 



GTliT, {pi}, {Ni}) - N Ru g R ^ 2 (T,p Rn0 , 



(3) 
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Here we also used that No — 2Ao 2 and = (1/2) M(V 
The Gibbs free energies of the slab and of the Ru02 bulk 
can be calculated by density-function theory, evaluating 
the total energies and the vibrational spectral As ap- 
parent from eq. 1(3}, r )(T,po 21 pco) depends only on the 
difference of the Gibbs free energies of the bulk and the 
slab results. In our previous work we could show that 
for Ru02 the vibrational energy and entropy contribu- 
tions to this difference of Gibbs free energies cancel to 
a large extend, and we will therefore replace the slab 
and bulk Gibbs free energies by the corresponding total 
energies. We note, however, that for other systems this 
procedure can cause a noticeable error. Furthermore, we 
note (details are discussed in the appendix) that some 
surface phases may be disordered which gives rise to con- 
figurational entropy. On the energy scale relevant to the 
present study, this configurational-entropy contribution 
is negligible (cf. the appendix). 



B. Chemical potentials 

The chemical potentials of O and CO, that enter eq. 
(|3*|). are determined by the condition of thermodynamic 
equilibrium with the surrounding gas phase reservoirs. 
Thus, their temperature and pressure dependence is 



Vo{T,po 2 ) 



and 



z^total 

h o 2 



fio 2 (T,p°) + k B Tln 



Po 2 

™0 



(4) 



(2iV Ru - N ) no(T,po2) - N C oHco(T,pco) 



Mco(I> co ) = E^ 1 + fi co {T,p a ) + k B Tln . 

(5) 

The temperature dependence of fio 2 (T,p°) and 
fLCo(T,p°) includes the contributions from vibra- 
tions and rotations of the molecules, as well as the ideal 
gas entropy at 1 atmosphere. It can be calculated, but 
in this paper we will simply use the experimental values 
from thermodynamic table o 9 i 21 , given in Table [I] 

The T = K values of the O2 and CO chemical poten- 
tials, Eo 2 tal and Eqq 1 , are the total energies of the iso- 
lated molecules (including zero point vibrations), and for 
these we will use our DFT results. We note that state-of- 
the-art DFT total energies of the O2 molecule are subject 
to a noticeable error. Again, for the present study this 
error is not crucial and can partially be circumvented 9 , 
however for other systems it maybe important to treat 
£totai differently (see , for example, the case of silver oxide 
discussed in Ref. |22|) . 

In the following we will present the resulting surface 
energies as a function of the chemical potentials, i.e. in 
(Mo> Mco)-space. Equations 10} and JSJ describe how 
these chemical potentials can be converted into pressure 
scales at any specific temperature. We will exemplify 
this below by showing pressure scales at T — 300 K and 
T = 600 K in order to elucidate the physical meaning 
behind the obtained results. 
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TABLE I: p,o 2 (T,p°) and p,co(T,p°) in the temperature 
range of interest to our study. The employed entropy and 
enthalpy changes are taken from the JANAF thermochemical 
tables at p° = 1 atm.— 



T Ao 2 (T,P°) 


Aco(T,P°) 


T ik) 2 (T,P ) 


Aco(T,P°) 


100 K 


-0.16 eV 


-0.14 eV 


600 K 


-1.22 eV 


-1.18 eV 


200 K 


-0.34 eV 


-0.33 eV 


700 K 


-1.46 eV 


-1.40 eV 


300 K 


-0.54 eV 


-0.53 eV 


800 K 


-1.70 eV 


-1.64 eV 


400 K 


-0.76 eV 


-0.73 eV 


900 K 


-1.98 eV 


-1.88 eV 


500 K 


-1.00 eV 


-0.95 eV 


1000 K 


-2.20 eV 


-2.12 eV 




FIG. 1: Sideview of the chosen supercell setup (Ru = light, 
large spheres, O = dark, medium spheres). Each slab consists 
of three 0-(RuO)-0 trilayers, and the atoms of one such tri- 
layer have been whitened in the lower slab. Consecutive slabs 
are separated by a vacuum region of ~ 11 A. 



using a (5 x 10 x 1) Monkhorst-Pack grid with 50 (15) k- 
points in the full (irreducible) part of the BZ. The energy 
cutoff for the plane wave representation in the interstitial 
region between the muffin tin spheres was -E™ f ax = 20 Ry 
for the wave functions and E™^ = 169 Ry for the poten- 
tial. With one notable exception this is exactly the same 
basis set as used in our preceding work on RuO2(110) 
in contact with a pure O environment^. The calculated 
short CO bondlength of 1.15 A (exp: 1.13 A) forced us 
to reduce the oxygen muffin-tin radius from 1.3bohr to 
1.1 bohr, so that we had to increase £™ f ax from 17 Ry to 
20 Ry to achieve the same high level of convergence as 
detailed before. 

With this basis set, the computed binding energies for 
the free 2 {E%» d = -6.07 eV), for the CO {E%g d = 
-11.42eV), and" for the C0 2 (E%fg = -17.62 eV) 
molecule are in very good agreement with previously 
reported DFT values, while still showing the known 
overbinding with respect to experimen1j2L22iSi. Note, 
that the above values contain zero-point vibrations that 
were estimated as 0.09 eV, 0.13 eV and 0.30 eV respec- 
tively. By comparison with test calculations performed 
at -E™ f ax = 24 Ry, we conclude that the absolute binding 
energies of O and CO at the RuO2(110) surface are con- 
verged within 0.15eV/atom, which translates into sur- 
face free energy variations below 10 me V/A 2 . For the 
differences between the computed 7(T, {pi}) of different 
surface phases, which are the quantities actually enter- 
ing into the construction of the phase diagram presented 
here, the error cancellation is even better. As a con- 
sequence, the numerical uncertainty of results discussed 
below is ±5meV/A 2 (with respect to the basis set and 
the supercell approach). While this error bar does not af- 
fect any of the physical conclusions drawn, we note that it 
does not include the errors introduced by the more basic 
deficiency of density-functional theory, namely the ap- 
proximate nature of the employed exchange-correlation 
functional, the effect of which will be discussed in sec- 
tion III C. 



C. DFT computations 

The total energies entering eq. J3J are obtained by 
DFT calculations using the full-potential linear aug- 
mented plane wave method (FP-LAPW) 23-24 - 25 together 
with the generalized gradient approximation (GGA) for 
the exchange-correlation functional 2 ^. The RuO2(110) 
surface is modeled in a supercell geometry, employing a 
symmetric slab consisting of three rutile 0-(RuO)-0 tri- 
layers as shown in Fig. ^ All atomic positions within 
the outermost trilayer were fully relaxed. A vacuum re- 
gion of « 11 A ensures the decoupling of the surfaces of 
consecutive slabs as described previously^ 

The FP-LAPW basis set parameters are: -R^t =1-8 
bohr, i?MT =1-1 bohr, i?£ T =1.0 bohr, wave function 
expansion inside the muffin tins up to l^ax = 12, poten- 
tial expansion up to = 4- For the RuO2(H0) slabs 
the (1 x 1) Brillouin zone (BZ) integration was performed 



III. RESULTS 
A. Stability range of R11O2 

As we are interested in the adsorption of reactants on 
a stable Ru0 2 substrate, we first analyze the stability 
range of rutile Ru0 2 bulk. In our previous workS we had 
already shown that in a pure oxygen environment, the 
/xo-variations may be restricted to a finite range: Below 
the so-called "O-poor limit" the oxide will decompose 
into Ru metal and oxygen, i.e., the oxide is only stable if 

A < ^ u u lk + 2Mo, (6) 

i.e., if 

Amo > \ [gZh 2 ~ SR U lk - £o° tal ] . (7) 
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Here A/z is defined as no -{l/2)E t ° t 2 ai . For T = K the 
right hand side equals half of the low temperature limit 
of the heat of formation, Hf(T = 0K,p — 0), for which 
our DFT calculations give H f (T = 0K,p = 0) = -3.4 
eV (the experimental result is —3.19 eV— ). Thus, despite 
the mentioned error in Eq* , for Hf there is a fortuitous 
error cancellation. We note in passing that for the bulk 
phases the T and p dependence of the Gibbs free energies 
is small. Thus, replacing the right side of eq. Q by 
(l/2)Hf(T = OK) is a good approximation. 

The "O-rich limit" refers to conditions where oxygen 
will condense, and this gives the other restriction: 

Afio < 0. (8) 

Combining the above two equations gives the range of 
A/io we will consider in the later discussion as 

±H f (T = OK, p = 0) < Afio(T,p 02 ) < 0. (9) 

Ru02 can also be destroyed (reduced) by carbon 
monoxide. In a pure CO environment the stability con- 
dition for the oxide is 

5ru0 2 + 2 Mco < 9 Z k + (10) 
This can be written as 

(11) 

Here A/ico is defined as nco — Eqq 1 . For the right hand 
side of eq. (JTTJ (at T = K) our DFT result is -1.5 eV 
(the experimental result is —1.33 eV—). 

Let us now consider that O2 and CO are both present 
in the gas phase. However, as mentioned earlier, we will 
consider the oxide (and its surface) in what we call a 
"constrained equilibrium" , i.e. we consider that R11O2 
is in thermodynamic equilibrium with O2 and with CO 
individually while the O2 and CO in the gas phase are 
not equilibrated. Combining eqs. Q and Ijllll defines a 
line in (A/io, A/ico )-space: 

A MO - A/ico = C, (12) 

with C = -0.2 eV (theory) and -0.26 eV (experiment 
Above the line, i.e. for A/ico > A/i — C, Ru02 decom- 
poses into Ru and below the line the oxide is stable. The 
situation close to the line will be discussed, also for other 
oxides, in a subsequent publication. 

Let us add a comment on the above treatment of the 
CO2 molecules. As mentioned above, CO2 is assumed to 
be not in thermodynamic equilibrium with any reservoir 
or any other species. A dynamical treatment of the CO2 
gas flow is outside the scope of this paper, and as a first 
estimate we therefore use the single-molecule internal en- 
ergy (total energy plus vibrations and rotations) in eq. 
(f 101) and in the following. 

Corresponding to the thus derived stability conditions 
we will only show our computed surface free energies in 




FIG. 2: Top view of the RuO2(110) surface explaining the lo- 
cation of the two prominent adsorption sites (bridge and cus), 
which continue the bulk-stacking sequence (central panel). 
Also shown are perspective views of the two most stable termi- 
nations involving only O (left panel) and involving O and CO 
(right panel). Ru = light, large spheres, O = dark, medium 
spheres, C = white, small spheres. Atoms lying in deeper 
layers have been whitened in the top view for clarity. 

the limited range -1.68 eV < A/i (T,p) < 0.0 eV and 
-2.0 eV < A/ico (T,p) < 0.0 eV, and mark the instabil- 
ity line, eq. 112JI . by a white dotted line. Finally, we also 
note that under very reducing conditions (low A/io, high 
A/ico) CO will transfrom into graphite, which then de- 
termines the C chemical potential and which e.g. occurs 
for A/ico = 0.0 eV for any A/iQ < — 1.2 eV. The corre- 
sponding region will be marked by extra hatching in our 
graphs, although it corresponds to gas phase conditions 
where Ru02 is already only metastable with respect to 
the aforedescribed CO-induced decomposition. 

B. Adsorption of O and CO on RuO 2 (110) 

In the rutile Ru02 bulk structure, every Ru atom has 
six oxygen neighbors and every oxygen atom three metal 
neighbors*^. The RuO2(110) surface structure is conve- 
niently explained on the basis of the mixed (RuO) plane 
termination exhibiting threefold coordinated O lattice 
oxygen atoms, cf. the middle panel in Fig. O The 
computed lattice constants for the rectangular (lxl) 
surface unit-cell are (6.43A x 3.12A), which compare 
well with the experimentally reported values of (6.35 A 
x 3.11 A)22i2i. The surface exhibits two distinct adsorp- 
tion sites, namely the so-called coordinatively unsatu- 
rated (cus) site atop of fivefold coordinated Ru atoms 
(Ru cus ), as well as the bridge site between two fourfold 
coordinated Ru atoms (Ru br ), cf. Fig. O The bulk stack- 
ing of the oxide would be continued by oxygen atoms first 
occupying all bridge sites (O br ), this way yielding the sto- 
ichiometric surface termination (upper left panel of Fig. 
[2Jl that also arises naturally when cutting the bulk be- 
tween two of the subsequent 0-(RuO)-0 trilayers shown 
in Fig. ^ Occupying also all cus sites (O cus ) then leads 
to the third possible (1 x 1) termination of a rutile (110) 
surface shown in the lower left panel of Fig. [21 9 

Adsorption of either O or CO is in principle also con- 
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TABLE II: Calculated binding energies in (1 x 1) phases on 
RuO2(110). The values are given with respect to the free 
O2 and CO molecule respectively, zero-point vibrations not 
included. 



Species 


Phase 


Binding Energy 


Qbr 


O br /- 


-2.44eV/atom 


Qbr 




-2.23 eV/atom 


Qbr 


o br /co cus 


-2.37 eV/ atom 


QCUS 




-0.99 eV/atom 


Q 3f 


o br /- 


-3.59 eV/atom 


co cus 


o br /co cus 


-1.38 eV/atom 


co cus 


co br /co cus 


-1.39 eV/atom 


co br 


co br /co cus 


-1.58 eV/atom 



ceivable at other sites on the surface. However, extensive 
calculations with O and CO in such low-symmetry sites 
showed significantly lower binding energies compared to 
adsorption at the cus and bridge sites, which is why we 
will restrict our discussion to these two prominent sites in 
the following. We also performed a large number of cal- 
culations in larger supercells up to (1x4) and (2 x 2) peri- 
odicity to assess the importance of lateral interactions be- 
tween the adsorbates. The obtained binding energies for 
both adsorbates varied by less than 150 meV—, reflecting 
rather weak lateral interactions beyond the direct first 
neighbor cus-br interaction (implicitly contained within 
our DFT-computcd total energies of (1 x 1) phases). Cor- 
respondingly only the latter phases will be considered in 
the following. 

Depending if the bridge and cus site within a (1 x 1) 
cell is either occupied by O or CO or empty, nine different 
adsorption geometries result. For these we introduce a 
short-hand notation indicating first the occupancy of the 
bridge and then of the cus site, e.g. O br /~ for O adsorp- 
tion at the bridge site, the cus site being empty. Having 
computed the total energies of all nine combinations, we 
only find the four depicted in perspective views in the 
side panels in Fig. [5] to be relevant in the context of a 
high pressure gas phase and will consequently dedicate 
ourselves to them for clarity. 

From our calculations we obtain a number of T = 
OK binding energies, which are given in Table [H] 
Good agreement for CO is obtained with other DFT 
calculations22i2&2!i2£. Also the value for O br agrees nicely 
with the value reported by Liu et alw&. On the other 
hand, all oxygen binding energies are ~ 0.8 eV stronger 
than those reported by Seitsonen and coworkers^MiS. 
We are unable to explain this discrepancy, but note that 
our O cus binding energy of about -1 eV (with respect to 
l/2i?Q 2 tal ) is consistent with the thermodesorption tem- 
perature of T = 300 — 550 K reported for this species^ 

From the numbers listed in Table [H] we see that O ad- 
sorption at bridge sites is in general significantly stronger 
than adsorption at the cus sites, which is comprehensible 
considering the on-top, onefold bonding geometry at the 
cus sites compared to the bridge, twofold bonding geom- 



Po 2 ( atm ) 

10 10 1 10 10 

H 1 h 600K 



10- 40 10- 30 10- 20 10 10 1 

H 1 1 1 1 300K 




Aju Q (eV) 

FIG. 3: Surface free energies of the three only O involving 
RuO2(110)-(l x 1) terminations (dashed lines) and of three 
O and CO involving surface geometries (solid lines). Corre- 
sponding geometries are explained in the left and right panels 
of Fig. 1, respectively. The dotted vertical lines indicate the 
allowed range of oxygen chemical potential, cf. eq. JUJ, while 
A/uco = 0.0 eV, corresponding to very CO-rich conditions. In 
the top x-axis, the dependence on A/io(T, po) has been cast 
into pressure scales at fixed temperatures of T = 300 K and 
T = 600 K. 

etry at the bridge sites, cf. Fig. |21 Interestingly, CO 
does not follow such a pronounced bond-order conserva- 
tion trend and exhibits only slight differences in bond 
strength at bridge and cus sites. As has already been 
pointed out by Seitsonen et at, this might be related 
to the specific CO br adsorption geometry, which in the 
here considered (lxl) phase is an asymmetric bridge 
with the CO much stronger bound to only one of the two 
neighboring Ru br atoms ^ 



C. Surface phase diagram 

Before proceeding to set up the surface phase diagram 
by evaluating eq. @ for any (/j.OjMco), we start by 
summarizing the previously reported effect of a just oxy- 
gen containing environment". In this special case, cor- 
responding to A/ico ~~ ► —00, we have to consider the 
stability of the aforedescribed three possible (lxl) ox- 
ide terminations -/-, O br /- and O br /O cus as a function 
of the oxygen chemical potential. From the computed 
surface free energies drawn as dashed lines in Fig. |3| 
we see that in the allowed range of A/^o delimited by 
the two vertical dotted lines, cf. eq. ©, the mixed 
(RuO) plane termination -/- shown in the topview in 
Fig. |21 is never stable, so that oxygen atoms will al- 
ways occupy at least all bridge sites, O br /-, as long as 
the surface is in equilibrium with the surrounding gas 
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CO br / CO cus 




Pa, (atm) 



J o (eV) 



FIG. 4: Surface free energies of the four stable geometries 
shown in Fig. 1, as a function of (A/io, A/xco). 



phase. Towards higher A/io oxygen atoms will addition- 
ally occupy all cus sites, leading to the polar O br /O cus 
termination, cf. Fig. |3J To give an impression of the 
corresponding pressures required to stabilize this latter 
termination, we have also included in Fig. pressure 
scales as second x-axis for T = 300 K and T = 600 K, 
the latter temperature corresponding to a typical an- 
nealing temperature frequently employed in experimen- 
tal studies on this system 14 i 15 i 16 i 17 i 19 i 33 i 34 i 35 . From these 
scales we deduce that the stoichiometric O br /- termina- 
tion, that was frequently prepared and studied in UHV 
studiea 14 i 15 i 16 i 17 i 19 i 33 i 34 i 35 , is primarily the end result of 
a high temperature anneal in UHV, whereas at atmo- 
spheric O2 pressures the surface exhibits additionally 
oxygen atoms at the cus sites. 

Extending our analysis now to the two-component 
(02,CO) environment, we fix the CO chemical poten- 
tial first to A/ico — 0.0 eV, corresponding to very CO- 
rich conditions. This allows us not only to conveniently 
draw the surface free energies of three most stable O and 
CO containing geometries into Fig. [3] (solid lines) , but 
gives also the low j(T, {pi})-limit for these phases. As 
explained in Section IIIA, in this CO-rich limit the C 
chemical potential is determined by the equilibrium with 
graphite below A/io = — 1.2 eV and with CO above, 
which is why the corresponding three lines in Fig. |3| 
exhibit a kink at this value. Even in the shown CO-rich 
limit, the CO br /- phase due to CO adsorption at only the 
energetically more favorable bridge sites, cf. Table [n] is 
barely more stable than the hitherto considered pure O- 
terminations. This changes, if CO is additionally present 
also at the cus sites, as we find the completely CO covered 
CO br /CO cus surface to have a low surface free energy at 
least at low O chemical potential. Towards higher A/zo 
the third geometry exhibiting O atoms at the bridge sites 
and CO at the cus sites, O br /CO cus , becomes finally even 
more stable, cf. Fig. |3J 

Allowing (in addition to A/io) variations of A/ico 
away from the CO-rich limit, the two-dimensional graph 
in Fig. is extended into the three-dimensional one 
shown in Fig. in which ~f(T, {pi}) is given as a function 
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FIG. 5: Surface phase diagram of stable and metastable struc- 
tures of RuO2(110) in (/io, A t co)-space. The additional axes 
at the top/right give the corresponding pressure scales at 
T = 300 K and 600 K. The white regions close to the bound- 
aries between different stable phases mark phase coexistence 
conditions, where at least 10% of the respective other phase 
is also present at the surface at T = 300 K. Regions marked 
by white hatching are particularly strongly affected by kinet- 
ics (see text). The white dotted line marks the stability limit 
of bulk Ru02 with respect to CO-induced decomposition: In 
the upper left part above this line, Ru02 is only metastable. 



of (A/io, A/xco)- This implies obviously, that the con- 
tents of Fig. [3]are contained as a special case in Fig. 0] in 
fact as the surface free energy dependence along the front 
x-axis of the latter figure. With decreasing CO chemical 
potential, i.e. lower CO content in the gas phase, the CO 
involving surface phases become increasingly less favor- 
able, whereas the surface free energies of the two pure 
O terminations do of course not depend on A/ico- It 
is also worth pointing out that according to eq. © the 
slope with respect to A/ico of the plane representing the 
CO br /CO cus phase in Fig. His twice that of the plane of 
the O br /CO cus phase, given that the prior structure con- 
tains double the amount of CO at the surface. Although 
A/xco could in principle be varied down to —00, we only 
show in Fig. 0]the range to — 3.0eV, as already at this 
value the pure O terminations have become most stable 
for any A/io, indicating that the CO content in the gas 
phase has become so low, that no CO can be stabilized 
at the surface anymore. 

Although Fig. ^ is instructive for understanding the 
extension of the surface free energy dependence from 
a one-component to a two-component environment, the 
three-dimensional nature of the plot does not allow an 
easy access to the really important information contained 
in it, i.e. which phase has the lowest j(T, {pi}) and is 
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thus most stable for a given gas phase characterized by 
A/iq and A^qq. This information is better obtained 
by only drawing the stability regions of the most stable 
phases, which corresponds to a projection of the lowest 
surface free energies in Fig. 0] onto the (A//o, A//co)- 
plane. The resulting surface phase diagram is drawn in 
Fig. [SJ At the lowest CO chemical potentials shown, 
so little CO is present in the gas phase, that we simply 
recover the previously discussed surface structure depen- 
dence on /iq in a pure O2 atmosphere. In other words, 
the stoichiometric O br /- phase is stable at low O chem- 
ical potential, while higher O2 pressures stabilize addi- 
tional oxygen at the surface, leading to the O br /O cus 
phase. Increasing the CO content in the gas phase, CO is 
first bound at the cus sites. This is easier at low O chem- 
ical potentials, where these sites are free, but is harder 
at higher A/j,o where also oxygen atoms compete for ad- 
sorption at the cus sites. As a result, towards higher A/io 
the O br / CO cus phase becomes only more stable than the 
O br /O cus phase at progressively higher CO chemical po- 
tentials, cf. Fig. [S] Finally, under very reducing condi- 
tions (high A/^cO: lo w A/io) CO is able to substitute O 
also at the bridge sites, yielding the completely CO cov- 
ered CO br /CO cus surface. Additionally shown in Fig. [5] 
is the line marking the instability of Ru02 with respect 
to CO-assisted decomposition at high A/ico (white dot- 
ted line). From the location of this line we see that the 
complete stability region of the CO br /CO cus phase (in- 
cluding the graphite-formation region hatched in black) 
and large parts of the O br /COcus phase correspond al- 
ready to metastable situations that will not prevail for 
long under realistic conditions. 

The surface phase diagram in Fig. [S] summarizes the 
key results of the present work. However, before pro- 
ceeding to discuss the physics contained in it, we like 
to discuss the accuracy of our theory. As stated in sec- 
tion IIC, the numerical uncertainty in the surface free 
energies due to the DFT basis set and the supercell ap- 
proach is about ±5 meV/A 2 . If we allow the four planes 
shown in Fig. |I]to shift by this value with respect to each 
other, the resulting intersections, i.e. the phase bound- 
aries, would shift within the white regions drawn in Fig. 
EJ We have added additional axes in Fig. [5] indicating 
the pressure scales at T — 300 K and 600 K. From these 
scales it becomes apparent, that at room temperature 
our computational uncertainty may very well correspond 
to even a few orders of magnitude in pressure. On the 
other hand, this is only so, because at T = 300 K these 
one or two orders of magnitude in pressure sample only 
a vanishingly small part of phase space. In general, we 
believe that Fig. illustrates quite nicely that our nu- 
merical inaccuracy does not affect the general structure 
of the obtained phase diagram at all - even if it would, 
it would just mean that we would have to increase our 
basis set. 

This error analysis does not yet include the more ba- 
sic deficiency of DFT, namely the approximate nature 
of the exchange-correlation functional. To this end, we 
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FIG. 6: Same as Fig. 4, but now computed using the LDA as 
exchange-correlation functional. 

have also computed the surface phase diagram not with a 
gradient-corrected functional, but with the local-density 
approximation (LDA)2&. The result (obtained by a com- 
plete geometry reoptimization of all structures using the 
LDA lattice constants) is shown in Fig. and is very 
similar to the GGA results (Fig. 03). It is worth point- 
ing out that the absolute surface free energies for the 
different phases obtained within the LDA are up to 20 
meV/A 2 higher than the GGA results. However, their 
differences, which ultimately determine the positions of 
the phase boundaries, differ only little. Although this 
LDA and GGA comparison does not provide a quanti- 
tative assessment of the remaining error that is due to 
the approximate treatment of exchange and correlation, 
the obtained small difference strongly suggests that the 
DFT-accuracy for this phase diagram is rather high. In 
line with our previous analysis^ we therefore conclude 
that the approximate nature of the exchange-correlation 
functional may affect the exact phase transition tempera- 
tures and pressure (with a typical uncertainty of ±100 K 
and 1-2 orders of magnitude in pressure), but does not af- 
fect the overall structure of the obtained phase diagram. 

D. Surface phase diagram, kinetics and catalysis 

The main outcome of the preceding section is the sur- 
face phase diagram of RuO2(110) in "constrained ther- 
modynamic equilibrium" with an environment formed of 
O2 and CO in the complete range of experimentally ac- 
cessible gas phase parameters. Kinetic effects due to 
catalytic CO oxidation at the surface may obviously 
cause deviations from this situation, and we will show 
in the following how the obtained surface phase diagram 
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may be analyzed to identify corresponding regions in 
(T, po2>-Pco) _s P ace where such kinetic effects may be- 
come crucial and/or a high catalytic activity may be ex- 
pected. 

In general, the surface will be close to our "con- 
strained thermodynamic equilibrium" with the reactant 
gas phase, as long as the on-surface CO2 formation is 
the rate-limiting step. Then, adsorption and desorption 
of the reactants can occur frequently on the time scale 
of the reaction event, allowing the chemical potentials 
of gas phase and adsorbed particles to equilibrate. In 
this context we note that the only proper definition of 
thermodynamic equilibrium is this equality of chemical 
potentials. This is not equivalent to saying that the ad- 
sorption and desorption rate of O and CO is equal. This 
would only be true, if no other channels for removing the 
particles from the surface existed. Instead, in catalysis 
we can have the situation, that desorption and reaction 
compete with adsorption, but this must not necessar- 
ily mean that the chemical potentials of the reactants 
couldn't still equilibrate with the ones in the gas phase 
(provided the reaction event rate is low). 

On the other hand, if the final reaction step is not 
rate-limiting and the bottleneck of the reaction would 
e.g. be given by the adsorption process(es), adsorbate(s) 
would be faster consumed by the final reaction than re- 
plenished by the slow adsorption from the gas phase. 
Under such conditions, the adsorbate concentration can 
be much lower or due to preferential adsorption much 
different to what is predicted by a thermodynamic the- 
ory. Still, in the present case the kinetic effects will even 
then be restricted to the immediate surface population at 
bridge and cus sites, provided one stays within the range 
of the Ru02 bulk stability discussed in section IIIA, i.e. 
below the dotted white line in Fig. [5] the creation of 
vacancies in the oxide bulk lattice is simply too costly, 
cf. the binding energy of 3f in Table |UJ A Mars -van- 
Krevelen type mechanism, where the oxidation reaction 
is noticeably dependent on oxygen diffusion through the 
bulk 37 , appears therefore rather unlikely for RuO2(H0) 
in a more oxidizing gas phase (A/io > A/ico)- Only if 
the steady-state reaction is run under gas phase condi- 
tions close to the dotted instability line in Fig. [5] would 
we expect this to become different, in which case also os- 
cillations between metallic and oxidic state might become 
perceivable^ 



E. Catalytically most active regions 

Analyzing the surface phase diagram in Fig. [3] along 
the lines described in the preceding section, the part of 
the O br /CO cus phase below the dotted instability line 
appears to be a likely region in (T,po,Pco)-space where 
catalysis might be most efficient: In contrast to all other 
phases, both reactants are then adsorbed at the oxide sur- 
face, indicating as dominant reaction mechanism CO cus 
+ O br -> C0 2 . We will show below that the C0 2 for- 



mation energy barrier for this process is noticeable (1.25 
eV, in good agreement with the result of Liu, Hu and 
AlaviSii). Furthermore, the energy to create vacancies in 
the O-bridge layer is rather high (2.37 eV, cf. TableU]}, 
so that we would expect the surface to be rather close to 
the "constrained equilibrium" situation under such gas 
phase conditions. 

The mechanism between O br and CO cus has al- 
ready been suggested on the basis of extensive UHV 
experimentsi^iS*i 7 *i2'i2iMiS ) where it had initially been 
characterized as a Mars-van-Krevelen type reaction, 
given that O br is adsorbed at the bridge sites which 
correspond to a continuation of the bulk oxide lattice 
stacking. As most of the experiments were not con- 
ducted at steady-state, but employed CO postdosage in 
UHV 15 i 16 i 17 i 18 i 34 i 35 , a further reduction of the surface was 
also reported, in which after O br had been reacted off, 
CO subsequently occupied also the bridge sitesi£'i2*i2* 3 £. 
While this is fully compatible with the higher binding 
energy of CO at the bridge sites, cf. Table ITU the phase 
diagram in Fig. ^reveals that under consideration of the 
environment, an extremely large CO/O2 partial pressure 
ratio (corresponding to an almost pure CO atmosphere) 
would be required to really be able to stabilize CO at 
the bridge sites in the CO br /CO cus phase, cf. the given 
pressure scales. Even then such a situation would not 
prevail for long under realistic conditions, as Ru02 is 
then already instable against CO-assisted decomposition 
into Ru metal. 

Hence, the initially identified part of the O br /CO cus 
phase is still a more likely candidate for a catalytically 
active region. The question remains if this is already the 
part in (T, p)-space where we expect the Ru02 catalysis 
to be most efficient? In this respect we note that a so- 
called stable phase in our thermodynamic approach is not 
stable on a microscopic scale, but represents an average 
over fluctuations of elementary processes such as dissoci- 
ation, adsorption, diffusion, association and desorption. 
As all these processes and their interplay^ are of cru- 
cial importance for catalysis, regions in (T, {pi})-space 
where such fluctuations are particularly pronounced are 
expected to be most important. This will be the case un- 
der gas phase conditions, which correspond to regions in 
the computed surface phase diagram close to boundaries 
between different stable phases: At finite temperatures, 
the transition from one phase to the other will not be 
abrupt in (A/io, A/zco)-space, but over a pressure range 
in which the other phase gradually becomes more pop- 
ulated. The resulting phase coexistence at the catalyst 
surface could then lead to a significantly enhanced dy- 
namics, in which even additional reaction mechanisms 
(in case of microscopically coexisting phases) or reaction 
fronts (in case of domain pattern formation) might be- 
come operational. 

Assuming a canonic distribution of the two compet- 
ing phases^ we have correspondingly estimated the region 
on both sides of the boundaries in which the respective 
other phase is present at least at a 10% concentration, 
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and marked it as white regions for the case of T = 300 K 
in Fig. [SJand Fig. Concomitantly, one of these bound- 
ary regions (the one between the O br /CO cus and the 
QbiyQcus phases, henceforth termed cus-boundary), for 
which we would thus expect an enhanced dynamics, falls 
just into gas phase regions, for which the exceptionally 
high turnover rates over working Ru catalysts have been 
reported, i.e. for ambient pressures in an about equal 
partial pressure ratioii*i£. Along this cus-boundary both 
reactants compete for adsorption at the cus sites. Pro- 
vided the rather low lateral interactions between neigh- 
boring cus sites described in section IIIB a microscopic 
coexistence of both phases appears then likely, i.e. neigh- 
boring cus sites could be occupied independently by ei- 
ther O or CO. This would open up an additional reaction 
channel, CO cus + O cus — > CO2, which had already been 
observed as a dominant mechanism in UHV experiments 
on preoxidized RuO2(110), as long as O cus atoms where 
available on these surfacesiii^. 

Below we report a rather low barrier of 0.9 eV for this 
process, so that O cus atoms can react away rapidly un- 
der these conditions. If the filling of empty sites with CO 
will also be fast as is the case in the hatched region of 
the O br /O cus phase close to the cus-boundary, then the 
surface will not be close to the constrained thermody- 
namic equilibrium situation and the real O cus concentra- 
tion may be much lower than suggested in Figs. EJandEl 
From our argument we therefore expect high catalytic ac- 
tivity in this region, but also note that here coverage and 
structure (i.e. the very dynamic behavior) must be mod- 
eled by statistical mechanics. Although the equilibrium 
approach thus breaks down just in this most interesting 
region, we see that it still enables us to rationalize un- 
der which gas phase parameters highest activity is to be 
expected. 

While this already helps to embed the UHV measure- 
ments into the catalytic context, the phase diagram also 
allows to systematically analyze when and how the pres- 
sure gap may be bridged in corresponding experiments: 
Namely, when it is assured that one stays within one 
phase region or along one particular phase boundary, as- 
suming that then the same reaction mechanism will pre- 
vail. For example, at an about equal partial pressure 
ratio of O2 and CO a decrease of ambient gas phase 
pressures over many orders of magnitude will at room 
temperature only result in a wandering close to the cus- 
boundary in (A/iOj A^co)-space, cf. Fig. Cor- 
respondingly, we would conclude, that similar reaction 
rates should result, as has indeed recently been noticed 
by Wang et ali^ when comparing their UHV steady-state 
kinetic data with high pressure experiments of Zang and 
KischSa. 

However, from Fig. |SJit is also obvious that without the 
knowledge of a (T, {pi})-phase diagram, as e.g. provided 
by the present work, a naive bridging of the pressure 
gap by simply maintaining an arbitrary constant partial 
pressure ratio of the reactants may easily lead to cross- 
ings to other phase regions and in turn to uncomparable 



results. In this respect it is particularly important to no- 
tice that e.g. the cus-boundary does not exactly fall on 
the equal partial pressure ratio line. There is no reason 
it should, and already at T = 600 K the same afore- 
described decrease of ambient gas phase pressures would 
bring one rapidly away from the cus-boundary until at 
UHV pressures one even ends up deep inside the stability 
region of the O br /- phase, cf. Fig. |S] Consequently, the 
computed phase diagram leads us to predict that kinetic 
measurements at this temperature and UHV pressures 
will result in a very low catalytic activity, as under such 
conditions no CO can bind to the surface anymore. This 
reactivity would then, however, be contrasted by the ex- 
tremely high turnover numbers reported for working Ru 
catalysts close to ambient pressures at exactly the same 
temperatur e - 11 ! 12 , in other words a typical pressure gap 
situation. 



F. Reaction mechanisms 

In the course of the preceding section two likely oxida- 
tion mechanisms have been identified, in which CO cus re- 
acts with oxygen adsorbed at either a neighboring bridge 
or a neighboring cus site. As apparent from Fig. [5] both 
channels take place along high symmetry lines parallel to 
the surface: Along the [110] direction for reaction with 
O br (dashed line) and along the [001] direction for reac- 
tion with O cus (dotted line), respectively, so that we are 
in both cases left with a planar configuration and may ne- 
glect the degrees of freedom perpendicular to this plane 
formed by the corresponding in-surface and the out-of- 
surfacc [110] direction. Even assuming the underlying 
oxide substrate as rigid, we are thus facing at least a six- 
dimensional problem, involving as one possible choice the 
following set of variables: Lateral and vertical position of 
C cus , lateral and vertical position of the adsorbed O, as 
well as the CO bondlength and the polar angle of the 
CO molecular axis with respect to the surface. In order 
to shed more light on the importance of the two possible 
reaction mechanisms we proceed to determine their re- 
spective transition state (TS) using large (1 x 2) unit-cells 
to decouple the periodic images of the reacting species. 
Within these cells, we map out the potential energy sur- 
face (PES) along two reaction coordinates, namely the 
lateral positions of the C cus and the adsorbed oxygen 
within the given reaction plane, while minimizing the en- 
ergy at each point with respect to all remaining degrees 
of freedom (including full a relaxation of the underlying 
oxide lattice). 

Figure[7]shows the mapped PES for the reaction CO cus 
+ O br — > CO2, which according to our phase diagram is 
likely for a large range of gas phase parameters corre- 
sponding to the stability region of the O br /CO cus phase. 
The total lateral distance between a neighboring cus and 
bridge site is 3.21 A, cf. Fig. |2 and we find the TS at 
a rather large lateral C cus displacement of 1.44 A, while 
the O br atoms move only by 0.21 A away from their equi- 
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C cus position along [110] (A) 



FIG. 7: Potential energy surface for the reaction CO cus + 
O br -» C0 2 . The lateral positions of C cus and O br along 
the [110] direction connecting the cus site with the neigh- 
boring bridge site, cf. the dashed line in Fig. 1, have been 
constrained, fully relaxing all remaining degrees of freedom. 
Actually calculated points are indicated by white circles. The 
energy zero corresponds to the initial state at (0.00 A, 3.21 
A) and the transition state geometry is shown in the inset 
(only the atoms lying in the reaction plane itself are drawn 
as threedimensional spheres). 



librium position. As shown in the inset of Fig. [7| the 
Ru cus -C cus distance is therefore remarkably lengthened 
at the TS (2.40 A compared to 1.99 A at the initial state), 
whereas the Ru br -O br distance is only slightly stretched 
(1.98 A compared to 1.92 A). Looking at the binding 
energies compiled in Table [n] we can identify as one rea- 
son for this asymmetric behavior the significantly larger 
binding energy for O br compared to CO cus . Nevertheless 
we compute a rather low reaction barrier of 1.25 eV for 
this process, which compares favorably with the value of 
1.15 eV reported in an earlier pseudopotential study by 
Liu, Hu and Alavi employing a constrained minimiza- 
tion for the O br -CO cus bond length^. As also pointed 
out in that work, already this mechanism on RuO2(H0) 
should have a higher reactivity than the Ru(0001) sur- 
face, for which the CO oxidation barrier was computed 
to be 1.45 eV in a p(2 x 2) ceYlM 

Turning to the second mechanism, CO cus + O cus — > 
CO2, we notice that the TS location in Fig. [8] has 
strongly shifted compared to the one displayed in Fig. \7\ 
At the TS also the O cus atoms have now moved laterally 
by 0.73 A stretching their bondlcngth from 1.70 A at the 
initial state to 1.79 A, so that the Ru cus -C cus bondlength 
needs only to be increased to 2.08 A in comparison to 
the 2.40 A required for the prior reaction. This differ- 
ence can again be understood by looking at the binding 
energies listed in Table [HI from where the much lower 
bond strength of O cus becomes apparent, that should fa- 




C cus position along [001] (A) 



FIG. 8: Potential energy surface for the reaction CO cus + 
qcus _^ C q 2 The lateral p 0sitions f C cus and qcus along the 

[001] direction connecting two neighboring cus sites, cf. the 
dotted line in Fig. 1, have been constrained, fully relaxing all 
remaining degrees of freedom. The actually calculated points 
are indicated by white circles. The energy zero corresponds 
to the initial state at (0.00 A, 3.12 A) and the transition state 
geometry is shown in the inset (only the atoms lying in the 
reaction plane itself are drawn as threedimensional spheres). 



cilitate an easier displacement in comparison to the much 
stronger bound O br atoms. In line with the thus ex- 
pected higher reactivity of the O cus atoms, we indeed 
find a much lower barrier of only 0.89 eV for this mecha- 
nism. 

In view of the surface phase diagram shown in Figs. 
H and this CO cus + O cus -> C0 2 reaction is primar- 
ily expected under the phase coexistence conditions along 
the cus-boundary. Concomitantly, it is also gas phase pa- 
rameters that correspond to this cus-boundary, for which 
the exceptionally high turnover numbers over working 
Ru catalysts have been reportedii*^. Obtaining a par- 
ticularly low barrier for a reaction mechanism that may 
preferably be operational under just these conditions fits 
thus not only nicely to the observed high reactivity, but 
also underlines as specific example, that the enhanced dy- 
namics in phase coexistence regions might be crucial to 
understand the function of surfaces under realistic con- 
ditions. 



IV. SUMMARY 

In conclusion we have computed a first-principles 
based phase diagram of the lowest energy surface struc- 
tures of Ru0 2 (HO) in the "constrained equilibrium" with 
an environment formed of O2 and CO. Depending on the 
chemical potentials of both gas phase species a number of 
different surface phases is found, ranging from two differ- 
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ent exclusively O containing terminations to a completely 
CO covered surface. We showed how this surface phase 
diagram can be analyzed to identify gas phase conditions 
where kinetic effects due to the CO oxidation reaction 
may become crucial and/or a high catalytic activity can 
be expected. In view of the possibly enhanced dynamics, 
we emphasize the particular importance of phase coexis- 
tence regions close to boundaries in the computed surface 
phase diagram, and suggest that a reliable bridging of the 
pressure gap is possible, provided that one stays within 
one phase region or along one particular phase boundary. 

Concerning the specific application to the CO oxida- 
tion over RuO2(110) we showed that the gas phase pa- 
rameters, for which exceptionally high turnover numbers 
have experimentally been reported, correspond to phase 
coexistence conditions at the catalyst surface, in which 
the reaction CO cus + O cus — > C0 2 may be active in ad- 
dition to the hitherto favored CO cus + O br -» C0 2 . For 
the prior mechanism we calculate a significantly lower 
barrier of 0.89 eV compared to the 1.25 eV for the latter 
channel, which for this specific example underlines that 
phase coexistence conditions might indeed be relevant to 
understand the reactivity of this catalyst surface. 

Our analysis exemplifies that calculations of the kind 
presented in this work can be used to identify important 
reaction steps at any given (T, {pi}) in the gas phase and 
to already explain a number of experimental findings on 
the surface reactivity. Yet, a real microscopic description 
of catalysis can only be obtained by subsequent compu- 
tations of the kinetics of the manyfold of possible ele- 
mentary processes, as well as simulations addressing the 
statistical interplay among them. To obtain meaningful 
results from the latter type of studies the consideration 
of all relevant atomistic steps is crucial, the identification 
of which will be greatly facilitated by the prerequisite of 
knowing the various stable surface phases in equilibrium 
with the given gas phase. 



Appendix 

Consider a system with N surface sites and a small 
number of n defect or adsorbate sites (n <C N). Then, 
the configurational entropy 5' config is given by 



N\ nl 



(13) 



If we define Ajitc as the surface area per site, the config- 
urational entropy per surface area is 
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(14) 



For N, n 3> 1 we can apply the Stirling formula which 
gives 
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The expression in the square brackets varies between 
(for (n/N) = 0) and 1.2 (for (n/N)=10%). Thus, the 
configurational entropy 



T5 .co„fi g k T 

< 1.2- 



Ait 



(16) 



For RuO2(110), considering two sites (bridge and cus) 
per (1 x 1) unit-cell, implies that A s { te = 10.03 A 2 . Cor- 
respondingly, we deduce a configurational contribution 
to the surface free energy of less than 5 meV/A 2 for any 
T < 1000 K. This is negligible everywhere in the phase 
diagram apart from the phase coexistence regions. 
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